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^\ . The magnetorotational instability (MRI) is thought to play a key role in the formation of stars 

' and black holes by sustaining the turbulence in hydrodynamically stable Keplerian accretion discs. 

I In previous experiments the MRI was observed in a liquid metal Taylor-Couette flow at moderate 

. Reynolds numbers by applying a helical magnetic field. The observation of this helical MRI (HMRI) 

< I ' was interfered with a significant Ekman pumping driven by solid end-caps that confined the insta- 

PL(' bility only to a part of the Taylor-Couette cell. This paper describes the observation of the HMRI 

, in an improved Taylor-Couette setup with the Ekman pumping significantly reduced by using split 

end-caps. The HMRI, which now spreads over the whole height of the cell, appears much sharper 
' and in better agreement with numerical predictions. By analyzing various parameter dependencies 

we conclude that the observed HMRI represents a self-sustained global instability rather than a 
noise-sustained convective one. 
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^ \ Cosmic magnetic fields are produced, on a wide range of spatial scales, by the homogeneous dynamo effect in 
^ • electrically conducting fluids [1]. The Earth's magnetic field is generated in its outer core, most likely by spiral flow 
structures which result from the combined action of Coriolis forces and thermal and/or compositional buoyancy 
The magnetic field of the Sun is generated in the convection zone and in the tachocline, very likely under the influence 
of differential rotation and helical turbulence 3] . A similar mechanism is probably at the root of large scale galactic 
magnetic fields Q , while inter-galactic magnetic fields are thought to be a product of the so-called fiuctuation dynamo 
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, Magnetic fields play an active role in cosmic structure formation. It was in 1991 that Balbus and Hawley lg| 
highlighted the importance of the magnetorotational instability (MRI, or Velikhov-Chandrasekhar instability [3, \^) 



for the angular momentum transport in accretion disks around stars and black holes. In MRI, the externally applied 
■ magnetic field serves only as a trigger for the instability that actually taps into the rotational energy of the flow. This 
is quite in contrast to another magnetic instability in which prevailing currents in the fluid can become unstable by 
. themselves. This latter, so-called Tayler instability 9] is held responsible for parts of the dynamo mechanism in stars 
' [l3|, and for some observed helical structures in jets and outflows [ll|. 
. ^ , Significant theoretical and computational progress has been made in the understanding of dynamo action and 
magnetic instabilities. Besides this, the last decade has also seen considerable experimental activities to investigate 
\ those effects in the liquid metal laboratory [l^, At the end of 1999, the hydromagnetic dynamo effect was 

. - - I observed nearly simultaneously at the large-scale liquid sodium facilities in Riga [l3| and Karslruhe [l^. In 2006, 
the VKS-experiment in Cadarache has shown not only self-excitation [l^ but also impressive magnetic field reversals 
[l3], although the use of soft-iron impellers with high relative permeability makes the interpretation of these results 
cumbersome. 

As for the so-called "standard MRI" (SMRI), with only an axial magnetic field applied, the authors of [3 claim to 
have observed it in a spherical Couette flow, but their background state was already fully turbulent, thereby defeating 
the original idea that the MRI would destabilize an otherwise stable flow. 

There is a simple reason for the fiow being turbulent before the SMRI sets in. The point is that the azimuthal 
magnetic field component of the MRI mode can only be produced from the axial field by induction effects, which are 
proportional to the magnetic Reynolds number (i?m) of the fiow. Rm, in turn, is proportional to the hydrodynamic 
Reynolds number according to Rm — PmRe, where the magnetic Prandtl number Pm — v jX is the ratio of viscosity 
V to magnetic diffusivity A — l/iJio<^- For liquid metals Pm is typically in the range 10~^...10~^. Therefore, in 
order to achieve Rm ~ 1, we need Re ^ 10^. ..10^, and wall-constrained flows with such high Re are known to be 
turbulent, whatever the linear stability analysis might tell. Apparently in contrast with this common wisdom, a 
group in Princeton has recently shown [l§] that a fiow profile in a short Taylor-Couette cell can be kept stable until 
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Re ^ 10^ if an appropriate configuration of split end-caps is used. However, the mechanisms that ultimately leads 
to flow stabilization in this configuration are not well understood. At least, the optimum rotation ratios of the split 
rings do not agree with numerical predictions, a discrepancy that might be explained by some wobbling of the inner 
cylinder [2^. 

Having learned that Rm 1 (and hence Re 10^) is needed just to produce the azimuthal magnetic field component 
of the MRI mode from the applied axial field, one might ask why not substitute this induction process by simply 
externally applying an azimuthal magnetic field as well? This idea was pursued in where it was shown that 

this "helical MRI" (HMRI), as we now call it, is indeed possible at far smaller Reynolds numbers and magnetic field 
amplitudes than SMRI. 

In [2l| it was also shown that HMRI and SMRI are continuously connected. In Fig. 1 of [2l[ the critical Rm was 
plotted against the rotation ratio ^ := font/ fin of outer to inner cylinder of a Taylor-Couette set-up. The extremely 
steep increase of this curve at the Rayleigh line, which occurs for a purely axial magnetic field, is just smeared out 
when an azimuthal field is added. Although HMRI, in the limit of small Re, is a weakly destabilized inertial oscillation 
[2^, there is a continuous and monotonic transition to SMRI when Re and the magnetic field strength are increased 
simultaneously. 

The relation of HMRI, which appears as a travelling wave, and SMRI is currently the subject of intense discussions 
in the literature [U, [2^, [2^, [l^, [28ll29l . [30l . [sil [s^ l , the roots of which trace back to an early dispute between Knobloch 
[33I [3^ and Hawley and Balbus [35l|. 

A remarkable property of HMRI for small Pm (which has been coined "inductionless MRI" ) , was clearly worked 
out in . It is the apparent paradox that a magnetic field is able to trigger an instability although the total energy 
dissipation of the system is larger than without this field. This is not so surprising though when seen in the context 
of other dissipation induced instabilities which are quite common in many areas of physics (36j . 

Another, and not completely resolved issue concerns the relevance of HMRI for astrophysical flows. On first glance, 
HMRI seems very attractive as it could extend the range of applicability of MRI into those regions of accretion disks 
which are characterized by a small Pm. This might apply to the "dead zones" of protoplanetary disks [s^] as well as 
to the outer parts of accretion disks around black holes [3g| . 

However, before discussing this point in more detail, it has to be checked whether HMRI works at all for Keplerian 
rotation profiles r2(r) ~ r^^/^ or not. In the inductionless limit, i.e. for Pm — 0, and using a local WKB analysis 
in the small-gap approximation, the answer to this question is "no" [23|. In contrast to this the solution of the 
corresponding global eigenvalue equation gives an affirmative answer, as long as at least the outer or the inner radial 
boundary is electrically conducting [s^. This is certainly good news for the working of HMRI in the outer, cold parts 
of accretion disks since their inner parts have indeed a higher conductivity. 

Unfortunately, even this is not the end of the story. Both methods, the WKB method and the global eigenvalue 
equation, deal with waves with a single wavenumber and frequency. However, since HMRI appears in the form of a 
travelling wave, one has to be careful with the interpretation of the positive growth rate of a single monochromatic 
wave. The crucial point here is that those monochromatic waves are typically not able to fulfill the axial boundary 
conditions at the ends of the considered region. To fulfill them, one has to consider wave packets. Only wave packets 
with vanishing group velocity will remain in the finite length system. Typically, the onset of this absolute instability, 
characterized by a zero growth rate and a zero group velocity, is harder to achieve than the convective instability 
of a monochromatic wave with zero growth rate. A comprehensive analysis of the relation of convective, absolute, 
and global instability can be found in the survey papers 40, 41]. In the context of MHD, this distinction played 
an important role, e.g., in the Riga dynamo experiment (42l |. A detailed analysis of the relation of convective and 
absolute instability for HMRI can be found in . From the extrapolation of the results of this paper it seems that 
Keplerian rotation profiles are indeed absolutely HMRI-unstable, but a final solution to this puzzle is still elusive. 

Notwithstanding these ongoing debates, the dramatic decrease of the critical rotation rates and magnetic field 
intensities for HMRI, as compared with SMRI, made this new type of MRI particularly attractive for experimental 
studies. 

First experimental evidence for HMRI was obtained in 2006 at the liquid metal facility PROMISE (Potsdam 
i?Ossendorf Magnetic /nS'tability iSxperiment) which is basically a Taylor-Couette (TC) cell made of concentric 
rotating copper walls, filled with GalnSn (a eutectic which is liquid at room temperatures). In (43| it was shown that 
the MRI travelling wave appears only in the predicted finite window of the magnetic field intensity, with a frequency of 
the travelling wave that was also in good accordance with numerical simulations. More detailed results were published 
in [4i,|4|. 

However, a slight difficulty with these early experiments was the fact that the MRI wave did not travel along the 
entire height of the Taylor-Couette cell but ceased to exist at some position typically close to mid- height. By modifying 
slightly the electric boundary conditions in radial direction, and analyzing the resulting change of the axial velocity 
[46|, it was possible to identify this sink of the travelling wave with the position of the radial jet that originates from 
the Ekman pumping at the upper and lower lids of the TC cell [43|. The influence of this Ekman pumping on the 
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flow structure was extensively discussed in a number of papers [H, [S^, . 

A related critical point is connected with electric currents which are induced in the lower end-cap that was made 
of copper in the first PROMISE experiments (now called PROMISE 1). In the worst case these currents are able to 
drive a significant Dean flow which possibly could modify the original Taylor-Coucttc flow proflle [27|, l48l|. 

In order to overcome both problems we have changed the axial boundary conditions, flrst by replacing the copper 
of the lower lid by insulating material, and second by splitting the (upper and lower) lids into two rings (rotating with 
the inner and outer cylinder, respectively) in order to minimize the Ekman pumping. Henceforth, this new version of 
the experiment will be called PROMISE 2 in order to discriminate it from the previous PROMISE 1 experiment. 

In this paper we will show that the minimization of the Ekman pumping makes the radial jet disappear almost 
completely. Consequently, the MRI wave can travel throughout the entire TC cell. We will investigate in detail four 
parameter dependencies of the MRI wave showing a much improved agreement with the 2D numerical predictions. 
Further comparisons with the results of a ID eigenvalue code for the convective and absolute instability will provide 
strong evidence that the observed wave represents a global instability and not only a noise-triggered convective 
instability as argued in [H, |3l| . 



II. FROM PROMISE 1 TO PROMISE 2 



Apart from the modified end-caps, the PROMISE 2 set-up is more or less identical to that of PROMISE 1 which was 
described in detail in [4^|4^. Its main part is a cylindrical vessel (cp. Fig.[T]) made of two concentric copper cylinders. 
The inner one is 10 mm thick extending in radius from 22 to 32 mm; the outer wall is 15 mm thick, extending from 80 
to 95 mm in r. This vessel is filled with the eutectic alloy Ga^^In^^-^Sn^^ '^ whose physical properties at 25°C are as 
follows: density p = 6.36 x 10"^ kg/m^, kinematic viscosity v = 3.40 x 10^^ m^/s, electrical conductivity a = 3.27 x 10^ 
(n m)~^. The magnetic Prandtl number is then Pm = ^^av = 1.40 x 10~^. 

The copper vessel is closed by a plastic bottom and is fixed, via an aluminum spacer, on a precision turntable; the 
outer wall of the vessel thus serves as the outer cylinder of the TC cell. The inner cylinder of the TC flow is flxed 
to an upper turntable, and is immersed into the liquid metal from above. It has a thickness of 4 mm, extending in 
radius from 36 to 40 mm. Between this immersed cylinder and the inner wall of the containment vessel there is thus 
a gap of 4 mm, also filled with GalnSn. The actual TC cell then extends in radial direction over a cylindrical gap of 
width d = Tout — '''in = 40 mm from the outer surface of the inner cylinder at rin = 40 mm to the inner surface of 
the outer wall of the containment vessel at rout = 80 mm, and in axial direction over the liquid metal height of mm 
< z < 400 mm. 

While in PROMISE 1 the upper end-cap was a plastic lid flxed to the outer frame this end-cap is now split into 
two concentric plastic rings (see Fig. [5]). The inner ring rotates with the inner copper cylinder, the outer one rotates 
with the outer cylinder. By this means, the opposing Ekman pumping effects at the rings compensate each other to 
a large extent. The choice of the splitting position at r = 56 mm was motivated by the simulation in (27j showing 
that a splitting at 0.4 d would minimize the global Ekman pumping for the case with the two rings attached to the 
cylinders (using only two rings significantly simplifies the construction when compared to a many-rings case). The 
lower end-cap is split in the same manner, so that we have a symmetric configuration, with respect to both the 
rotation rate and to the electrical conductivities. This is a significant improvement compared with PROMISE 1, in 
which the copper bottom was simply part of the vessel which rotated with the outer cylinder. 

The magnetic field configuration is identical to that of PROMISE 1, with axial magnetic fields of order 10 mT 
being produced by the current Icon in a double-layer coil with 78 windings, and an azimuthal flcld of the same order 
being generated by a current /rod through a central water-cooled copper rod of radius 15 mm. 

As in PROMISE 1, the measuring instrumentation consists exclusively of two high- focus ultrasonic transducers 
(from Signal Processing SA) with a working frequency of 4 MHz which are fixed into the outer plastic ring, 12 mm 
away from the outer copper wall, flush mounted at the interface to the GalnSn. Since this outer ring is rotating, it 
is necessary to transfer the signal into the laboratory frame (in contrast to PROMISE 1 in which the upper lid was 
fixed to the frame) . This is accomplished by the use of a slip ring contact which is situated below the TC vessel (not 
shown in Fig. [T|). 

The advantage of the Ultrasound Dopplcr Vclocimctry (UDV) is that it provides full profiles of the axial velocity 
Vz along the beam-lines parallel to the axis of rotation. The spatial resolution in axial direction is chosen 0.685 mm, 
the time resolution is 1.84 sec. The width of the beam (over which Vz is averaged) is approximately 8 mm, according 
to the diameter of the ultrasonic transducers. 

To get reliable signals from along the 400 mm height of the TC cell one has to ensure a good acoustic contact between 
the transducer and the GalnSn which is often compromised by non-perfect wetting, sticking oxides or bubbles. A 
good signal is all the more important since the velocity perturbations to be measured are typically in the order of 0.1 
mm/s, which is at the lower edge of the applicability of the UDV technique. 
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This fact has to be kept in mind when comparing, i.e., rms values of the velocity perturbation from different runs. 
For example, the identification of the transition between stable and unstable regimes rely strongly on the measured 
velocity distribution. However, this velocity signal is typically affected by noise which, in turn, depends on the acoustic 
contact. 

Typically, the duration of the experimental runs was 1900 sec, after a waiting time of one hour or even longer in 
cases when long transients were to be expected. 

III. RESULTS 

In this section, we will present the main exp erimental results and compare them with numerical results based on 
a ID eigenvalue solver (for details consult [UH^]) and on a time-stepping 2D solver that was already used in [27j . 
Usually, for both solvers perfectly conducting boundary conditions were applied at the inner and outer radius, which 
is certainly not completely correct since the conductivity of copper is approximately 13 times higher than that of 
GalnSn, but not infinite. For some parameters, at least, it has been checked that the inclusion of the finite-thickness 
copper walls in the ID eigenvalue code does not change significantly the results based on perfectly conducting walls. 

In particular we will study the dependence of the HMRI on the four parameters that govern the problem. Given 
the inner and outer radii of the TC cell, and rout (with a gap width d — r^ut — Tin), the rotation rate of inner and 
outer cylinder, /in and /out, the axial magnetic field and and the azimuthal magnetic field at the inner cylinder, 
B^p{r — Tin), these governing parameters are the Reynolds number 

Re := 2nfi^d/i^ , (1) 

the ratio of rotation rates 

:= /out//in , (2) 

the Hartmann number 

Ha := B,(ri„da/pi.)i/2 , (3) 

and the ratio of azimuthal to axial magnetic field 

/?:=S^(r = n„)/B, . (4) 

With the actual geometry and the material properties of GalnSn given, we can derive the following useful relations 
between the experimental parameters and the dimensionless parameters: 

Re = 29586 • /i„/Hz, (5) 
Ha = 0.158 • /eoii/A, (6) 

P = 0.0491 • /,od//coii. (7) 

Before entering into the details of the parameter dependencies of the MRI, we present the results in the magnetic 
field free case. 

A. Results without magnetic fields 

The fiow structure in the non-magnetic case is an important reference point for the later analysis of the MRI. 
Starting with /i = 0, we expect a Taylor vortex flow (TVF) whose amplitude decreases for increasing ^. Beyond the 
Rayleigh point /iRayi := ''hi/^out = 0.25, the TVF is supposed to disappear completely in accordance with Rayleigh's 
criterion. However, in any real TC-cell with finite height this is far from trivial due to the existence of a meridional 
flow driven by the Ekman pumping at the upper and lower lids. Approximately at mid-height of the TC cell the two 
opposed vertical flows meet each other and produce a radial jet [47[. 

The difference of the flows in the PROMISE 1 and PROMISE 2 conflguration, respectively, is documented in Fig. [31 
What is actually shown is the velocity component measured along the ultrasonic beam-lines. This velocity component 
is supposed to represent the total axial velocity Vz{z,t), although some slight influence of the strongly dominant 
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component cannot be ruled out in case the direction of the UDV transducers is not perfectly vertical. On the r.h.s. 
of each panel the time-averaged value va^iz) is given in dependence on the axial position z. 

Figure [3] shows, at /in — 0.1 Hz, the results for PROMISE 2, (a), (b), and (c) corresponding to ^ — 0, ^ — 0.1, and 
fi = 0.27, respectively. As expected, we observe a TVF which is getting weaker for increasing /i and which disappears 
completely at /i = 0.27. For comparison, Fig. [3ji shows one result for PROMISE 1, at /in — 0.06 Hz and /i — 0.27. 
Here we observe a very steep sign change of Vz{z) approximately at mid- height which necessarily must be connected 
with a radial jet at this point, a feature that is nearly absent in PROMISE 2 (Fig. Evidently the global effect 
of the Ekman pumping is drastically reduced by the splitting of the end-caps. In the following we will see that this 
absence of the radial jet has important consequences for the propagation of the MRI wave. 

B. Variation of Re 

As with many other instabilities in hydrodynamics, MRI is expected to set in at a certain critical Reynolds number 
-Rccrit- In order to study this transition in detail, we have fixed the rotation ratio to /x = 0.27 and the current through 
the coil to /coil = 76 A (i.e. Ha = 12). The current through the central rod has been set to either /rod = 4000 A 
(i.e. (3 — 2.59) or /rod = 7000 A (i.e. (3 = 4.53). Throughout the paper (if not otherwise stated) we have chosen the 
rotation to be clockwise, when viewed from above, and Bz and /rod directed downward, where the latter induces 
in the clockwise direction. 

To start with, we show in Fig. [3] the axial velocity signals, again the time averaged vaz{z) on the r.h.s, but here 
and henceforth the deviation Vz{z,t) :— Vz{z,t) — vaz{z) of the total velocity from the time averaged value on the 
l.h.s. The panels (a) and (b) exhibit the individual signals from the UDV sensors 1 and 2, respectively, (c) shows their 
average, and (d) shows their difference. Actually, in this difference signal one observes a slight m — 1 mode with a 
dominant frequency equal to the rotation rate of the outer cylinder, which results just from some slight asymmetries 
of the facility (e.g. the current leads to the central rods). This is in contrast to PROMISE 1 where we have observed, 
in some parameter regions jisl . |45| . [46| , an m = 1 mode with a frequency between the rotation rate of inner and outer 
cylinder whose origin is not yet clarified. 

For /rod = 7000 A and /coii ~ 76 A, the 12 panels of Fig. [5]show Vz{z, t) for increasing /in (i.e. Re). From here on, 
the l.h.s. of each panel exhibits the average over the two sensors (corresponding to panel (c) in Fig. |4]). Evidently, 
the travelling wave instability sets in approximately between /in = 0.04. ..0.05 Hz (i.e. Re = 1183. ..1479) and remains 
until the highest value /in = 0.24 Hz (i.e. Re = 6830) that we have studied in the experiment. 

It is also instructive to see what happens when we change the direction of Bz by changing the direction of the 
current in the coil. It is well-known that the direction of propagation, i.e. whether the wave pattern drifts in the +z 
or — z direction, depends on whether the product flBzB^p is negative or positive, respectively. Figure [H documents 
this effect for the parameter choice /in = 0.1 Hz, /x = 0.27, /rod = 7000 A, and IcoH = ±76 A. Note that the shape of 
the (a) upward and (b) downward travelling wave is now quite symmetric, in contrast to PROMISE 1 in which the 
position of the radial jet made the pattern more asymmetric (see, e.g. Fig. 4 in [46|). 

When speaking about upward and downward travelling wave it is also necessary to say that these two types do not 
appear exclusively. Actually, the upward travelling wave, when reflected at the upper boundary, produces a slight 
downward travelling wave, and vice versa. To analyze this point we show in Fig. [7] the 2D Fourier transform from 
the z,t-space to the fc, /-space of the signal Vz{z,t) for /in = 0.05 Hz. Not surprisingly, the experimental data are 
more smeared out than the numerical predictions. Nevertheless, the position of the maximum in the k, / space is 
approximately the same in both cases. The maximum of the 2D FFT is mainly localized in the second and fourth 
quadrant of the k, /-space which indicates an upward travelling wave. However, the slight signal in the first and third 
quadrant indicates also the existence of a downward travelling wave. 

By applying a Hilbert transform to the original data, it is possible to disentangle the dominant upward travelling 
wave and the minor downward travelling wave which results from a reflection of the upward travelling wave at the 
upper end-cap. This is done, now again for /in = 0.1 Hz, fj, = 0.27, in Fig. [H Inspired by [5(|, we will later use the 
ratio of the maxima of the downward to the upward travelling wave as an additional indicator for the onset of the 
absolute instability. 

In Fig. [5] we restrict ourselves to the downward travelling wave and compare its shape with simulated results for the 
particular cases /in = 0.04 Hz and /in — 0.05 Hz. We observe a quite convincing agreement of experimental data and 
numerical simulation, even the amazing slow decay of a transient wave at /in = 0.04 Hz is recovered by the numerics. 

A more quantitative comparison of various experimentally and numerically obtained features of the MRI wave, in 
dependence on /in, is made in Fig. \W\ This figure contains four parts: Fig. [TUb shows the squared rms of Vz{z,t), 
calculated between 250 mm< z <350 mm (this particular range has been chosen since the velocity in this interval is 
comparably high and the UDV signal is usually not much affected by noise as it is for smaller values of z). Figure 
[TUb shows the squared ratio of the maximum amplitude (taken over all z) of the downward to that of the upward 
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travelling wave. This quantity had been used successfully in other problems to distinguish between absolute and 
convective instability |50]. Figures [TUb and [TUtl show the normalized wave frequency /wave //in and the normalized 
phase velocity of the wave Cwavc/(/innn), respectively. The same dependencies, but now in more detail for 0.02 Hz 
< /in < 0.1 Hz are shown in Fig. [TTl 

Basically, the agreement between experiment and numerical results from the 2D code is quite convincing, in par- 
ticular the onset of the instability approximately at /in ~ 0.4 Hz for /rod = 7000 A. However, there is one significant 
deviation: while the numerics indicate a disappearance of the instability approximately at /in = 0.11 Hz, the exper- 
iment shows a continuation of the instability far beyond this point. The reason of this difference is not yet clear. 
What can be said, at least, is that the chosen parameter n = 0.27 is very close to the maximum value of /x beyond 
which HMRI ceases to exist. Here, the instability is extremely sensitive, and any slight deviation of the experimental 
conditions from the numerical assumptions can play a significant role. 

For /rod — 4000 A, the amplitude of the velocity is very small in good accordance with the numerics which predicts 
stability for all values of Re. 

As noted in the introduction, the importance of MRI lies in the fact that it destabilizes rotational flows with 
outward increasing angular momentum which would otherwise be hydrodynamically stable, according to Rayleigh's 
criterion. Doing experiments on MRI it is necessary to validate that the angular momentum of the base flow is indeed 
increasing outward. This is all the more important since a peculiarity of HMRI at small Reynolds and Hartmann 
numbers is that it works only for values of fi which are slightly above the Raylcigh line. Since, in a finite system, 
some deviation from the ideal TC-flow profile is unavoidable, in particular close to the end-caps, one has to be careful 
that the observed instability is indeed a magnetically triggered one and not an instability of a flow which is (at least 
in some parts) hydrodynamically unstable. 

Unfortunately, the UDV technique allows us only to measure axial velocity distributions in detail. Some efforts 
have been made to measure by means of an electric potential probe inserted into the liquid metal, but neither the 
necessary accuracy nor the desired spatial resolution could be achieved by this method. For that reason we have to 
rely on the numerical simulations exclusively. Given the good correspondence of numerics and experiment, as seen 
so far in the Vz data (at least for not too high values of Re), we believe that the numerics gives a more or less true 
picture of the real physical process. 

In the following we will consider the particular case n — 0.27, Icon = 76 A, /rod = 7000 A with two different values 
of /in = 0.04 Hz (corresponding to Re = 1183), and /in = 0.05 Hz (corresponding to Re = 1479). For these two cases 
we plot in Fig.[T2]the numerical snapshots of the angular velocity il(r, z), of the specific angular momentum r^fl{r, z), 
and of the radial derivative of this specific angular momentum. 

While the plots for ll(r, z) (Fig. [T2k ) are not very instructive and seem to indicate only a slight difference between 
the two cases, the plots for the angular momentum and its radial derivative show considerably more structure. For 
/in = 0.04 Hz we observe a rather regular outward increase of the angular momentum, with the only exception close 
to the splitting positions of the upper and lower end-caps. This is not surprising since, with the inner ring rotating 
with the (fast) inner cylinder and the outer ring rotating with the (slow) outer cylinder, there must be a certain region 
where the angular momentum decreases, while in the rest of the fluid it increases outward. 

The crucial point here is the transition from global stability to global instability when going from /in = 0.04 Hz 
to /in = 0.05 Hz. For the latter we observe a wave that travels throughout the total volume of the TC cell. Clearly 
visible is the "fingering" of the angular momentum which is typical for MRI fFll. 

C. Variation of /i 

Without any magnetic field, the case fi = corresponds to the classical TC flow with its well-known critical 
Reynolds number Recin = 68.2 (for our particular case Hn/fout = 0.5) (s^ . With a purely axial or a purely azimuthal 
field applied this critical value would increase. If both field are applied, the character of the instability changes 
from a steady TVF to a travelling wave, although with a very small frequency. If is then increased, the wave 
frequency increases also. HMRI means that the travelling wave continues into the region beyond the Rayleigh point 
MRayi '■= '''?n/''out = 0.25 where any purely hydrodynamic instability would disappear according to Rayleigh's criterion. 

For the PROMISE 1 experiment the behaviour of the wave with increasing /i was documented in Fig. 5 of [45*1 . It 
was characterized by a rather broad range of instability extending to values as large as fi — 0.35. Here we present 
the corresponding results for the PROMISE 2 experiment. Figure [T51 shows 12 panels for < /i < 0.35. We have 
fixed /in = 0.1 Hz, /coii = 76 (i.e. Ha = 12). The current through the central rod has been set to /rod = 7000 A (i.e. 
P = 4.53). It is clearly visible that the travelling wave continues beyond the Rayleigh line, though not very far. At 
/i = 0.32 it has disappeared completely. 

The comparison of the upward travelling wave for /i = 0.26 and /j, = 0.27 is illustrated in Fig. [TH Figures [T51 and 
[TBI show again the squared rms of the velocity, the squared ratio of the upward and downward travelling wave, the 
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wave frequency and the wave speed. In addition to the results of a time-stepping 2D solver we have also added the 
critical fi values for the convective and the absolute instability as they result from a ID eigenvalue solver [33 |. 

Besides the good agreement of experimental and numerical results from the 2D code, it is most remarkable that the 
transition point from the 2D code corresponds quite accurately to the absolute instability from the ID simulation, 
while the ID results for the convective instability lie at much too high /i values. This is strong evidence that the 
observed instability is a global one and not only a noise-triggered convective one as claimed in [31j . 

D. Variation of /3 

For liquid metals with Pm — 10^^.. .10"^, the critical Reynolds number for SMRI (with a purely axial magnetic 
field) is in the order of several million. The key point of HMRI is that by adding an azimuthal magnetic field, this 
number drops drastically to approximately 10'^. For PROMISE 1, the transition to instability for increasing f] turned 
out to be not very sharp, as documented in [4^ . In the following we will investigate the corresponding transition for 
PROMISE 2. 

Figure [T71 shows 12 panels for < /rod < 8200 A. Evidently, there is no indication of any instability for /j-od < 3000 
A. The transition to MRI occurs approximately between /rod = 4000 and /rod — 5000. The comparison with the 
numerical simulation for /rod = 4500 A and /rod — 7000 A is shown in Fig. [T51 

Figure [191 shows again the quantitative analysis for two different lines corresponding to fi = 0.26 and fi = 0.27. The 
first noticeable fact is that for /i = 0.26 we obtain, quite similar to the case of the /i line discussed above, a good 
agreement of the ID numerical results for the absolute instability and the 2D simulation, while the ID code results 
for the convective instability would point to much smaller values of the critical (3. Quite interesting, however, is the 
fact that for fi = 0.27 the agreement of absolute instability (at /rod ^ 8200 A) and global instability (/rod ^ 7000 
A) becomes worse. Evidently, the value /x = 0.27 is (at least for the given Ha number) very close to the upper 
limit of the HMRI. In this case, the threshold of MRI will be extremely sensitive to any slight modification of the 
base V^{r) profile. This is a typical case in which we indeed observe a global instability for the experimental axial 
boundary conditions, but the hypothetical use of ideal Taylor-Couette boundary conditions would make the fiow 
stable again. This interesting case has been observed in the numerical simulations of Liu ([1^ for PROMISE 1, [sij 
for PROMISE 2). However, as we see here the interpretation of the author is not correct: what is observed for real 
boundary conditions is not a noise- induced convective instability (this would start already at f3 — 2.3) but still a 
global instability, although on the basis of a flow that is slightly modified by the lids which make the fiow a bit more 
prone to MRI than the ideal TC-flow would be. 

E. Variation of Ha 

MRI is a weak field instability characterized not only by a lower critical Ha but also by an upper one. For the 
PROMISE 1 experiment it had been shown in [43| that the travelling wave appears approximately at a coil current 
30 A and disappears again at 110 A. 

Here we discuss the corresponding behaviour for PROMISE 2. We choose /i„ = 0.06 Hz, /i = 0.27, and study 
in detail the cases /rod = 7000 A and /rod = 8000 A. Figures [20l [211 and [22| document again the MRI dependence 
on /coil (i.e. Ha), the comparison with numerics, and the quantitative analysis, respectively. While the transition 
between stable and unstable regime for large Ha numbers (around 15) shows a satisfactory agreement of experiment 
and numerics, this agreement becomes worse for the transition at low Ha numbers. Actually, the observed MRI starts 
already at lower Ha (around 4. ..5) in contrast to the theory which predicts it to start between 8 and 9. This effect 
needs further work to be understood. 

Similar as for the /3-line (at fi = 0.27) in the former subsection, we have to note that the absolute instability would 
only start slightly beyond the chosen rod-current of 8000 A. We obtain, nevertheless a global instability which points 
to a slight modification of the base fiow. 

IV. CONCLUSIONS 

In this paper we have documented and analyzed four types of parameter variations in order to characterize in detail 
the HMRI in a TC-flow with split end-caps (PROMISE 2). Due to the strongly reduced Ekman pumping, resulting 
from this apparently minor modification, we have observed much sharper transitions from the stable fiow regime to 
the MRI regime than in the former PROMISE 1 experiment. 
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In general, the experimentally observed thresholds turned out to be in good correspondence with the thresholds 
resulting from 2D simulations for the global instability. In addition to that, we have shown that the results of a ID 
eigenvalue solver for the absolute instability are in good agreement with the global instability, at least as long as we 
are not too close to some upper limit values of /i (approximately 0.27 for the considered values of Re, Ha and /?) for 
the occurrence of HMRI. 

In contrast to that, the convective instability resulting from the ID code would generally imply a larger range of 
the instability than observed in the experiment. From this it seems justified to conclude that the observed instability 
is indeed a global one which corresponds, in a wide range of parameters, to the absolute instability and not to a 
noise-driven convective one. 

Coming closer to those parameters (in particular /i) for which helical MRI ceases to exist (at least for the achievable 
values of Ha and /3) we observe increasing deviations between global and absolute instability, which very likely trace 
back to the slightly changed velocity profile of the base flow in the finite TC-cell. 

What is left for future work is a careful comparison of the z-dependence of the wave amplitude with the corresponding 
results of the 2D code and of the ID eigenvalue code (in which this z-dependence appears as an imaginary part of 
the wavenumber). Comparisons for a few parameters seem to indicate a fairly good agreement, but a quantitative 
comparison for all parameter variations requires much more work. 

A "promising" perspective of PROMISE would be to extend the range of parameters Re and Ha in order to 
investigate the continuous transition between the well-studied region of HMRI and the experimentally much more 
demanding region of SMRI. One way to accomplish this would be to replace the working fluid GalnSn by sodium. 
Another interesting extension of the work would be to study various combinations of MRI and Tayler instability, 
which would require to direct electric currents through the liquid metal. The most interesting idea, however, namely 
to realize in experiment a self-sustaining nonlinear dynamo process in an MRI-destabilized (quasi)-Keplerian shear 
flow jssj , will probably remain a dream, at least for liquid metal flows. Some appropriate plasma experiments [s^l 
may be better suited for such an ambitious project. 
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FIG. 1: Sketch of the PROMISE 2 experiment. The gap width of the Taylor-Couette cell is 40 mm, the height is 400 mm. 




FIG. 2: Details of the PROMISE 2 experiment at the top. The inner ring rotates with the inner cylinder, the outer ring rotates 
with the outer cylinder. The same applies to the end-cap at the bottom. 
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FIG. 3: Measured axial velocity Vz{z,t) (l.h.s.), and time averaged vaz{z) (r.h.s.), in the field-free case for PROMISE 2, at 
fin = 0.1 Hz: (a) n = 0, (h) n = 0.1, (c) ju = 0.27. (d) The same for PROMISE 1 at /in = 0.06 Hz and n = 0.27. 
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FIG. 4; Measured axial velocity for fin = 0.1 Hz, /out = 0.027 Hz, I,od = 7000 A, and Icou = 76 A. (a) at UDV sensor 1, (b) at 
UDV sensor 2, (c) average of sensor 1 and sensor 2, (d) difference of sensor 1 and sensor 2. The frequency seen in (d) is just 
the rotation rate of the outer cylinder. 
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FIG. 5: Vz{z,t) and 110^(2;), both averaged over the two UDV sensors, in dependence on /in (i.e. Re) for n = 0.27, /rod = 7000 
A, and Zcoii = 76 A. 
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FIG. 6: Vz{z,t) and vaz{z), both averaged over the two UDV sensors, for /in = 0.1 Hz, p = 0.27, /rod = 7000 A: (a) 7coii = 76 
A, (b) /coil = —76 A. The change of the current direction in the coil makes the MRI wave change its direction from (a) upward 
to (b) downward. 
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FIG. 7: Power spectral density of the 2D FFT of Vz{z,t) for /in = 0.05 Hz, /i = 0.27, /rod = 
panel: experimental results, Lower panel: results of the simulation with a 2D solver. 
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FIG. 8: vaz{z) and Hilbert transforms of Vz{z,t) for /in = 0.1 Hz, jj, = 0.27, /rod = 7000 A, Icoii = 76 A. (a) Dominant upward 
travelling component, (b) minor downward travelling component which results from a reflection of the upward travelling 
component at the upper end-cap. 
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FIG. 9: Comparison of experimental and simulated results for the upward travelling wave component for /j, = 0.27, /rod = 7000 
A, /coil = 76 A, /in = 0.04 Hz (upper two panels) and /in = 0.05 Hz (lower two panels). 
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FIG. 10: Squared rms of Vz{z,t), squared ratio of maximum amplitude of downward to upward travelling waves, normalized 
frequency and normalized wave speed for /x = 0.27, /rod = 4000 A and 7000 A, Icon = 76 A, in dependence on /in (i.e. Re). 




FIG. 11: Same as Fig. 1101 but zoomed to medium values of /i 
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FIG. 12: From top to bottom: Simulated angular velocity fi(r, z), specific angular momentum r'^fl{r,z), and radial derivative 
of r^Q{r, z) for n = 0.27, /rod = 7000 A, Icoii = 76 A. (a) Subcritical /in = 0.04 Hz. (b) Supercritical /in = 0.05 Hz. 
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FIG. 13: Vz{z,t) and vaz{z), both averaged over the two UDV sensors, in dependence on n for /in = 0.1 Hz, I^od = 7000 
and /coil = 76 A. 
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FIG. 14: Comparison of experimental and simulated results for the upward travelling wave component for /i 
^rod = 7000 A, Jcoii = 76 A, yti = 0.26 (upper two panels) and n = 0.27 (lower two panels). 



= 0.1 Hz, 



23 




FIG. 15: Squared rms of Vz{z,t), squared ratio of maximum amplitude of downward to upward travelling waves, normalized 
frequency and normalized wave speed for /in = 0.1 Hz, I^od = 4000 A and 7000 A, Icon = 76 A, in dependence in fi. 
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FIG. 16: Same as Fig. 1151 but zoomed to fi values close to the Rayleigh point /iRayi 
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FIG. 17: Vz{z,t) and vaz{z), both averaged over the two UDV sensors, in dependence on /? for /in = 0.06 Hz, = 0.26, and 
/coil = 76 A. 
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FIG. 18: Comparison of experimental and simulated results for the upward travelling wave components for /in = 0.06 Hz, 
H = 0.26, and Jcoii = 76 A, /rod ~ 4500 A (upper two panels) and /rod ~ 7000 A (lower two panels). 
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FIG. 19: Squared rms of Vz{z,t), squared ratio of maximum amplitude of downward to upward travelling waves, normalized 
frequency and normalized wave speed for /in = 0.06 Hz, /j, = 0.26 and 0.27, /coii = 76 A, in dependence on /rod (i-e. (3). 
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FIG. 20: in dependence on I^oU (i-e. JJo) for /in = 0.06 Hz, fx = 0.27, and /rod = 8000 A. 
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FIG. 21: Comparison of experimental and sinmlatcd results for the upward travelling wave components for /in = 0.06 Hz, 
H = 0.27, and hod = 8000 A, Zcoii = 70 A and Icon ~ HO A. 
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FIG. 22: Squared rms of Vz{z,t), squared ratio of maximum amplitude of downward to upward travelling waves, normalized 
frequency and normalized wave speed for /in = 0.06 Hz, /j, = 0.27, /rod = 7000 A and 8000 A, in depencence on Icon (i-e. Ha). 



